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Abstract 

We present our results of Monte-Carlo simulations of polarized Compton X- 
rays from magnetic cataclysmic variables, with realistic density, temperature 
and velocity structures in the accretion flow. Our study has shown that the 
X-ray linear polarization may reach about 8% for systems with high accre- 
tion rates viewed at a high viewing inclination angle. This value is roughly 
twice the maximum value obtained by previous studies which assumed a 
cold, static emission region with a uniform density. We also investigate the 
X-ray polarization properties of ultra-compact double-degenerate binaries 
for the unipolar-inductor and direct-impact accretor models. Our study has 
shown negligible X-ray polarization for the unipolar-induction model. How- 
ever, the direct-impact accretor model may give X-ray polarization levels 
similar to that predicted for the magnetic cataclysmic variables. 



1.1 Introduction 

Magnetic cataclysmic variables (mCVs) and Ultra-compact double degener- 
ate binaries (UCDs) are potential X-ray polarization sources. mCVs contain 
a magnetic white dwarf accreting material from a low- mass, Roche- lobe fill- 
ing companion star. There are two major types: (i) the AM Herculis binaries 
(AM Hers, also known as polars) and (ii) the intermediate polars (IPs) (see 
Hi)). In AM Hers, the white-dwarf magnetic field {B ~ lO"^ - 10^ G) is 
strong enough to lock the whole system into synchronous rotation. It also 
prohibits the formation of an accretion disk, and the accretion flow is chan- 
nelled by the magnetic field into the magnetic polar regions of the white 
dwarf. The white dwarf in an IP has a weaker magnetic field {B 10^ G). 
The white-dwarf magnetosphere truncates the inner part of the accretion 
disk, and the material flow is channelled by the magnetic field from the 
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Fig. 1.1. An illustration of the emission region in mCVs where the polarized X-rays 
are generated. The accretion column has a circular cross-section, with a radius Tc. 
Compton scatterings occur within a height H from the white-dwarf surface. 

inner disk rim to the white-dwarf surface. For both AM Her and IP, the 
supersonic accretion flow becomes subsonic abruptly near the white dwarf 
surface, thereby forming an accretion shock. The shock heats up the accret- 
ing matter to keV temperature, and X-rays and optical/IR radiation are 
emitted from the shock-heated gas as it cools and settles onto the white- 
dwarf surface (see e.g. (EgI)). 

A UCD consists of two white dwarfs revolving around each other in a 
very tight orbit (which has a linear size of Jupiter). There are debates on 
what mechanisms generate the X -ray s in UCDs. The leading models are 
the unipolar-inductor (UI) model (|l9l ) and the direct-impact accretor (DIA) 
model In the UI model, the UCD orbital dynamics are determined 
jointly by magnetic interaction and gravitational radiation. Magnetic in- 
duction similar to that in Jupiter and lo sets up an electric current circuit 
across the binary. The electrical dissipation at the foot-points of magnetic 
field lines connecting the two white dwarfs heats up the white-dwarf atmo- 
sphere. Because of the convergent field configuration, small hot spots are 
formed at the surface of the magnetic white dwarf and emit X-rays. In 
the DIA model, the X-rays are accretion powered as in conventional binary 
X-ray sources. The close proximity between the two white dwarfs creates 
a strong tidal interaction and this prohibit the formation of an accretion 
disk. Mass transfer between the two stars is instead facilitated by a dense 
material stream. X-rays are emitted from a hot spot, where the accretion 
stream impacts onto the white-dwarf surface. The accretion configuration 
of the DIA model is somewhat similar to that of AM Hers. 
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The optical/IR emission from AM Hers is strongly polarized. The circular 
polarization can be as be as high as several ten percents (113) • The optical/IR 
polarization is generated by a thermal cyclotron process ([l|), where energetic 
electrons in the shock-heated plasma gyrate around the magnetic field. Op- 
tical/IR circular polarization has been observed in several IPs (see but 
it is much weaker than that of the AM Hers. The origin of the polarization 
is believed to be the same as that of AM Hers. In the UI model, UCDs are 
predicted to be strong sources of electron-cyclotron masers (jl4l : The 
masers are generated through a loss-cone or a shell instability (see (9|)) de- 
veloped in the electron population when the charged particles stream along 
the converging magnetic field lines near the polar regions of the magnetic 
white dwarf. The masers are characterized by a high brightness temperature 
and almost 100% circular polarization. 

Polarization in the low-energy emission of mCVs is produced by processes 
involving gyration of electron around a magnetic field. Their X-ray polar- 
ization is, however, unrelated to magnetic plasma processes. Instead it is 
due to scattering, where unpolarized thermal X-ray photons emitted from 
the shock-heated region are scattered by the electrons that precipitate onto 
the white dwarf. In this article we show the results of our calculations of X- 
ray polarization in mCVs and discuss briefly the astrophysical implications 
of our findings. We also show some preliminary results from our study of 
X-ray polarization in ultra-compact double degenerate binaries. 



1.2 Results and Discussions 

1.2.1 X-ray polarization in magnetic cataclysmic variables 

The geometrical setup in our calculations is shown in Figure II. 1[ The den- 
sity, temperature and velocity structure in the post-shock emission region 
is determined by the hydrodynamics model given in 
wise stated, thermal free- free emission is the dominant cooling process. We 
consider a Monte-Carlo approach to simulate the Compton scattering events 
and the transport of the polarized X-ray photons in the accretion column. 
A non-linear algorithm based on H; 0) is used to determine the photon 
mean-free path and the scattering probability. The scattering variables are 
determined following the prescriptions described in |l3 ). The polarizations 
are calculated using the Klein-Nishina cross section and the formulation for 
photon-electron scattering given in The polarized photons that escape 

from the accretion column are binned according to their energy-momentum 
and summed to give the spectral polarization at specific intervals of viewing 
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Fig. 1.2. Linear polarization Pl of X-rays as a function of the viewing inclination 
angle 6 for accretion onto magnetic white dwarfs with masses of 1 Mq. Filled 
triangles correspond to the cases with a specific accretion rate m = 1 g cm~^ s~^; 
and open circles correspond to the case with 10 times higher specific accretion rates. 
The Thomson scattering optical depths across the accretion shock are r = 0.04 and 
0.3 respectively. The H/vc ratios are 3, 5 and 10 (panels from top to bottom). 



inclination angles. The formulation, computational algorithms and numeri- 
cal simulation procedures are presented in detail in (j8|). 

Figures 11.21 and 1 1.31 show the results of two example simulations. For the 
case with a 1.0-Mq white dwarf, the maximum value of the linear polariza- 
tion Pl is about 1—2%. For the case with a O.S-M© white dwarf, the max- 
imum Pl may reach about 8% for viewing angle ~ 90°. The polarization 
increases slightly with 6 for low specific accretion rates (m ~ 1 g cm~^ s~^). 
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Fig. 1.3. Same as Figure fL^ for white dwarfs with masses of 0.5 Mq. The Thomson 
scattering optical depths across the accretion shock are r = 0.1 and 1.0 for the cases 
with low and high specific accretion rates respectively. 



However, for sufficiently high m (~ 10 g cm^^ s^^), the polarization could 
increase substantially at large 6. 

Figure 11.41 shows a comparison of the polarization from a structured ac- 
cretion flow and a cold, static and uniform density emission region. For low 
accretion rates, the cold, static, uniform density emission region underpre- 
dicts the linear polarization, especially at large 6. For high accretion rates, 
it overpredicts the polarization for angles below about 70 — 80°. The main 
difference of the two models is that the structured flow model always has 
a highly dense base region, where most of the scatterings occur, but the 
static model has a uniform density throughout the scattering region. The 
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Fig. 1.4. Comparison between the X-ray polarizations predicted by a static, cold 
emission region with a uniform density (solid lines) and by an emission region with 
structures determined by the hydrodynamic model of (|l8r ) (dotted lines). The 
parameters of the models are the same as those of the case in the middle panel 
of Figure 11751 Triangles correspond to cases with m = 1 g cm~^ s~^; and circles 
correspond to cases with m = 10 g cm~^ s~^. 
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Fig. 1.5. Comparison between X-ray polarizations from mCVs with accretion flows 
dominated by thermal free-free cooling and by cyclotron cooling. The white-dwarf 
mass is 0.5 Mq, and the specific accretion rate is 10 g cm~^ s^^. The ratios 
of efficiencies of cyclotron to thermal free-free cooling at the shock es — (open 
circles) and 10 (filled triangles). 

underprediction of the polarization by the uniform model at low rh is due 
to the lack of a dense base layer which gives a substantial scattering optical 
depth. The overprediction of the polarization by the uniform model at high 
m is caused by a uniformly high scattering optical depth at all heights above 
the white-dwarf surface, in contrast to the density drop off with height in 
the stratified accretion flow. 

The high density at the base and rapid density drop off with height are 
also the reasons why the angle dependence of the polarization is insensitive 
to the accretion-column aspect ratio in the cases with stratified accretion 
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Fig. 1.6. Predicted X-ray linear polarization in the mCV GK Per as a function of 
the white-dwarf spin phase (j). In the simulations the white-dwarf mass is set to be 
0.63 M0 flC^, and the specific accretion rate is assumed to be 10 g cm~^ s~^. 
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Fig. 1.7. Linear polarization as a function of viewing inclination angle for a UI UCD 
model. In the model, the bulk Lorentz parameter of the fast streaming electrons 
r = 10^, the effective Thomson scattering depth across the column t — 10~^. The 
H/rc ratio is 100. 

flows. The situation is different for a static, uniform-density scattering re- 
gion (see (0)), where the effective scattering optical depth depends strongly 
on the viewing inclination and the aspect ratio of the accretion column. 
Note the polarization is insensitive to the cyclotron cooling process for the 
same reason, as scatterings occur mainly in the dense base not the less dense 
region immediately below the shock (Figure [L5]) . Figure [L6l shows the pre- 
dicted polarization from a system with the same parameters as those derived 
for the IP GK Per. 



1.2.2 X-ray polarization in ultra-compact double degenerates 

The accretion geometry of the DIA UCD is similar to that of the AM Hers, 
except that an accretion shock might not be formed. In terms of Compton 
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scattering, the model set up would be approximately the same as that of 
the case with a static, uniform density accretion column. Although the 
gas has a bulk motion in the accretion column, the relatively slow speed 
implies that Compton recoil is the dominant effect in the scattering event. 
As a first approximation the consequence is not much different from that 
for the case with cold electrons. The polarization is therefore of the order of 
1 — 4% (see Thus, as a rough estimate, a polarization of a few percent 
would be expected from the DIA UCD model. The UI UCD does not have 
a dense accretion column. Instead there is a stream of relativistic electrons. 
The dominant effect is Doppler shift instead of recoil. Moreover, almost 
all events are head-on as seen by the relativistic electrons. Figure UCD 
shows our simulation for a model UI UCD. As the scattering optical depth 
is small, there is no substantial polarization. X-ray polarization can clearly 
distinguish the DIA and the UI model for UCD, despite the fact that the 
two models could have very similar X-ray spectral and timing properties. 
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